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Disclaimer 

Although this program was tested by its author prior to submission, 
no warranty, express or implied, is made by the author, 1620 USERS 
Group, or IBM as to the accuracy and functioning of the program and 
related program material and no responsibility is assumed by the 
author, 1620 USERS Group, or IBM in connection therewith. 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 



1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter or pencil, do not use ink) 

Program No. Date 

Program Name: 

1. Does the abstract adequately describe what the program is and what Yes No 

it does? 
Commen t 

2. Does the program do what the abstract says ? Yes No. 

Commen t 

3. Is the Description clear, understandable, and adequate? Yes No. 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No. 

Comment 

Are the Sense Switch options adequately described (if applicable)? Yes No. 

Are the mnemonic labels identified or sufficiently understandable? Yes No. 

Commen t 

5. Does the source program compile satisfactorily (if applicable)? Yes No. 

Comment 

6. Does the object program run satisfactorily? Yes No. 

C omment 

7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No. 

C omment 

8. Does the Program Meet the minimal standards of the 1620 Users 

Group? Yes No. 

C omment 

9. Were all necessary parts of the program received? Yes No. 

Comment 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 



Mr. Richard L. Pratt Company 
Data Corporation 

7500 Old Xenia Pike Address 
Dayton, Ohio 45432 



User Group Code. 



THIS REVIEW FORM IS PART OF THE 1620 USER GROUP ORGANIZATION'S PROGRAM 
REVIEW AND EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVITED 
TO PARTICIPATE IN THIS EVALUATION. 
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RELOCATION OFFSETS 



by 

P. 0. Roberts 
Assistant Professor of Civil Engineering 
M.l.T. Civil Engineering Systems Laboratory 



Purpose 

Frequently during the course of designing a highway or railroad, 
situations arise in which the design engineer would like to change the 
alignment that has been chosen. This can occur because of a misjudgment 
of the quantities incurred or because of a change in the design situ- 
ation, such as the continued growth of a subdivision or opposition of 
the public to the destruction of a park or historic building* Frequently, 
the change involved is small enough that the original survey information 
is still valid and usable. 

Whatever the reason for the change, the work involved in recom- 
putation is formidable. New geometry must be calculated. Equivalent 
stations on the new centerline must be calculated and offset distance 
from the original survey line, or base line, to the new centerline must 
be established. These offsets are then used for correcting the cross 
section sheets. The program is designed to relieve the burden of 
computation which alignment changes place on the engineer. See Figure 1. 

Description 

The program is accomplished in two parts. Part One computes the 
geometry of the baseline and the new centerline, given the state plane 
coordinates of the P.k's and one parameter defining each curve. Part Two 
works with the table of geometry prepared by Part One and computes the 
centerline stations and offsets for each baseline station. 
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Program Operation - Part One 



The geometry computation is performed as follows: 
J. The point identification, state plane coordinates of the 
P. I., and either the radius, the degree of curvature, the 
tangent length, or the external distance is given for each 
P.I. on the baseline. 

2. The cards containing this information are read in and stored 
in memory. The geometry computations are performed and punched 
out P. I. by P. i. 

3* The program then cycles back to the beginning and reads the 
same information for each P.I. on the centerline. 

k. As the output for each P.I. is computed and output, a geometry 
table is built up in memory. The P.C., P.T„, and P.I. stations, 
the azimuth of each line, and the deflection angle at each turning 
point are output for use by the engineer. 

Program Operation - Part Two 

The second program operates with the geometry table left in memory 
by program one and with a starting baseline station number and increment 
of advance. 

1. The baseline station number is used to compute the coordinates 
of that point on the baseline and the azimuth of the cross section 
located at that point. 

2. The line originating at the point where the cross section 
crosses the baseline and with an azimuth normal to the baseline is 
intersected with the centerline. The station of intersection offset 
and skew are determined. This information is output for each incre- 
ment as designated on the input. 
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PART 2. 



- Features 



1. The station number of the origin is placed in the radius field 
on the input card. 

2. Only one of the curve parameters need be specified. The others 
are entered as zero. If more than one is entered, the first is used 

3. The termini of both baseline and center line alignments are in- 
dicated by placing 999999-99 in the radius field of the input card. 

k. Several baselines and/or centerlines can be handled contin- 
uously. Intermediate termini of each line are indicated by 
888888.88 in the radius field. 

5> After an intermediate terminus, the radius field of the 
following card is expected to contain the origin station number 
of the new line. If this number is zero, the terminus station 
of the previous line is assigned. 

6. The starting station and the increment of advance of the 
stations on the baseline are specified on the first input card 
to part two. The stationing is automatically increased by the 
amount of the increment until the stationing exceeds the value 
on the next card, or until the baseline numbers fail to match. 
When this happens, the values on the second card replace those in 
memory. 

Restrictions 

1. A total of 30 P.I. points can be handled. P.I. points include 
origins and termini as well as P. I.'s on both the centerl ine and 
base! ine. 

2. All values on the input cards must be entered as either a 
number or as zero* They may not be left blank. 
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3. Accuracy is determined by the number of significant digits in 
the input. Eight significant digits are accepted by the program. 

k. Each P.I. must be numbered with a four digit ID number. The 
units and tens position should contain the point number (which is 
not used by the program). The hundreds and thousands digits must 
contain the baseline or centerl ine number. These numbers must 
match those on the part two input cards. 

5- Situations in which the baseline and centerl ine are going in 
opposite directions should be avoided. The program will not es- 
tablish an intersection when the azimuths of the two lines in- 
volved are more than 90 degrees apart. 

Operating Instructions 

1. PARITY switch on STOP. 

2. I/O switch on STOP. 

3. Flow switch on PROGRAM. 

k. Program Switches 

Switch 1 on for typed output 

Margins - 14 and 95 

Tabs- 22, 34, 46, 58, 70, 82 

5. Press - RESET. 

6. Ready read hopper with Part One Program followed by geometry 
input cards. 

7. Press - LOAD. 

8. Ready punch hopper with blanks. 

9. Press - PUNCH START. 

10. When program has read all cards and punched output, it will 
stop with a 48 in the operation register. 

11. Press - RESET. 

12. Ready read hopper with Part Two Program followed by the 
Baseline Station input cards. 

13- Press - LOAD. 

14. Program will stop on a 37 read instruction in the operation 
register after punching output. 
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RELOCATION OFFSETS OUTPUT DATA 



POINT 

101 
101 
102 
102 
103 
10$ 
104 
104 
201 
201 
202 
202 
301 
301 
302 
302 
303 
303 

4oi 

401 
402 
402 
403 
403 



GEOMETRY OUTPUT 
X/PC Y/PI 



500.00000 
00000000 
600.00000 
312.31042 
500.00000 
820.60377 
500.00000 
1420.6036 
550.00000 
2000.0000 
500.00000 
2304.1379 
600.00000 
1000.0000 
700.00000 
1168.6172 
600.00000 
1814.8315 
550.00000 
1814.8315 
400.00000 
2056.5387 
400.00000 
2732.8722 



100.00000 
.00000000 
500.00000 
412.31041 
900.00000 
820.60377 
1500.0000 
1420.6036 
1500.0000 
2000.0000 
1800.0000 
2304.1379 
100.00000 
1000.0000 
500.00000 
1412.3104 
900.00000 
1814.8315 
900.00000 
1814.8315 
1400.0000 
2336.8465 
1800.0000 
2732.8722 



R/PT 

.00000000 
.00000000 
359. 99993 
508.29335 
.00000000 
820.60377 
888868.88 
1420.6036 
2000.0000 
2000.0000 
999999.99 
2304.1379 
1000.0000 
1000.0000 
974.77287 
1646.2143 
888888.88 
1814.8315 
1814.8315 
1814.8315 
1909.8593 
2613.1801 
999999.99 
2732.8722 



D/AZ 

.00000000 
14.036246 
14.323947 
345.96376 
.00000000 
.00000000 



.00000000 
350.53768 



.00000000 
14.036246 
5.8778595 
345.96376 



.00000000 
343.30076 
3.0000000 
.00000000 



T/AMG 

.00000000 
.00000000 
99.999994 
•28.072491 
.00000000 
14.036246 



.00000000 
.00000000 



.00000000 
.00000000 
243.69325 
•28.072491 



.00000000 
.00000000 
280.30781 
16.699245 



S/CLENG 

.00000000 
.00000000 
12.310558 
195.98293 
.00000000 
.00000000 



.00000000 
.00000000 



.00000000 
.00000000 

30.000000 
477.59718 



.00000000 
.00000000 
20.460514 
556.64148 



SOURCE LANGUAGE LISTING 



OFFSETS OUTPUT 

CL STA OFF COS ANG 

1075.7464 97.01.4230 1 .0000000 

1 175.7490 96.988070 .99997324 

1275.9644 9U 109430 .99394231 

1394.9944 82.391260 .99991534 

1516.2619 88.469440 .99366932 

1618.4773 96.619620 .99959515 

1718.4814 97.014238 1.0000000 

1848.5383 48.193600 .99982753 

1897.7233 26.182357 .95782552 

2002 1264 -3.8179427 .95782552 

2106.3430 -33.139600 .96499239 

2209.2707 -57.446400 .97771940 

2311*0171 -76.144900 .98751011 

2411.9016 -89.402200 .99445141 

2419.3218 -140.43650 .99797536 

2429.3468 -139.77230 .99762776 

2439.3705 -139.05580 .99725273 

2449*3992 .-138. 28730 .99685002 

2459.4393 -137.46500 .99641932 

2469.4779 -136.59000 .99596115 

2479.5207 -135.66200 .99547526 

2489.5717 -134.68100 .99496141 

2499.6241 -133.64680- .99441993 

2509.6747 -I32.55CSO .993^5099 

2519.7393 -131.4ir90 .99325370 



BL 


BL STA 




100.00000 




200.00000 




300.00000 




400.00000 




500.00000 




600.00000 




700.00000 




800.00000 




900.00000 




1000.0000 




1 100.0000 




1200.0000 




1300.0000 




1400.0000 


2 


2010.0000 


2 


2020.0000 


2 


2030.0000 


2 


2040.0000 


2 


2050.0000 


2 


2060.0000 


2 


2070.0000 


2 


2080.0000 


2 


2090.0000 


2 


2100.0000 


2 


2110.0000 



/y 



O 



C RELOCATION OFFSETS / PAHT 1 

C P. 0. -ROP-ERTS / JULY* 1961 

C MIT CIVIL FNGINFFRINC SYSTry. LABORATORY 

C SW 1 ON FOR TYPED OUTPUT / KAH'l4t95 / TAtiS 22 . 34, <*6, 58 , 70 , 82 

DIMENSION ID ( 31 ) . V ( 30 ) . X < -0 V .R ( 30 ), PC ( 30 ) ,P I ( 30 ) .PT HO ) ,OFFD ( 3 ) 
DIMENSION D( 30 ) ♦ T f 30 ) »F ( 30 ) 
C . INITIALIZE 

100 l=n 
C READ ROUTINE 

DO 262 N=l»2 
ID(3l)=J 
DO 120 J»l »30 
1 = 1+1 

REAo,!ou)»xtn.Y(n»Rm»^i)'T<n.rm 

IF(R( 11-999999.98)1 ?0» 1^0.1 *f) 
120 CONTINUE 
140 I=I-J+1 

r . 

200 ISWlsl 
ISW3=1 
RI=R( H 
PHI)«RI- 

pet n«Ri 

PT ( I )=RI 

C=0.n 

T(I)=0. 

201 

GO TO 10 
206 IF( IS*1-2>248.210»?K> 
248 ISWW 

P 



ANG 10=0.0 
ALENG=0. 
GO TO 232 

210 ANGI = AZ-AZO 

212 IF(ANGI-3. 1415927) 211.211.214 

211 ANGI=ANGI+6. 2831854 
GO TO 212 

214 ANGI* ANGI -6 .2 83 18 54 
IF(ANGI)208. 209.209 

208 DS=-1 .0000000 
GO TO 205 

209 DS=1. 0000000 
205 ANG=ANGI/2. 

C FIND WHICH IS GIVEN 

I F ( R ( I ) )221, 216.221 

216 IF<D(I))221»217,219 

217 IF(T( I ) )221 .218.220 

218 IF ( E ( I ) ) 221 .222.223 
C IF D IS GIVEN 

219 R ( I ) =5729. 5780/D( I ) 
GO TO 221 

C IF T IS GIVEN 

220 R ( I ) =T ( I )*COS( ANr,)/SIN(ANr,)#ns 

GO TO 221 
C IF E IS GIVEN 

223 R ( l ) =E ( I )/(! ./COS( AMG)-1 . ) 

GO TO 221 
C IF NONE ARE GIVEN 



222 D«I)=0.0 

GO TO 215 
221 DU|=5729.5780/R( I ) 
215 T(n=Rtl )*tSIN<A.NG)/COS(ANG) )*DS 

PC( I )=PI( I.)-T( I ) 

ALENG*ANGI*R( I )*DS 

PTt I )=PCt ! )+ALENG 

EC I l»l 1«./CQS( ANG )-l • )*RC I ) 

ANGIO=ANGI*57. 295780 

232 A20UT» A2*57.295780 

I F( SENSE SWITCH 1)233*234 

233 PRINT » IDM ) *X ( 1 ) # Y I I > »R < I ) *D ( I ) »T ( I ) »E ( I ) 
PRINT* ID* I » *PCC I ) »P I ( I ) »PT ( I ) » AZOUT *ANGIO» ALENG 
GO TO 235 

234 PUNCH* I0( I ) »X( I ) ,Y( T ) »P( I ) »H( I ) »T ( I ) »F( T ) 
PUNCH. IDtf )»PCC I ) »P J( I ) »PTM J * AZOUT , A*IG!0»ALENG 

235 PHI+l)«PTtn-T(I)+C 
A20=AZ 

1*1 + 1 

240 ' IFtRI I )-888888«B7)?42»244.744 
242 CONTINUE 
GO TO 201 

244 PCin-PI (U 
PT( I )=PI ( I I 

fF(SENSE SWITCH 1)245*246 

245 PR I NT » I D ( I ) » X M ) » Y ( I ) » 3 ( I ) 
PRINT>IDfI)*PC(I)*PI(I)»PT(I ) 

GO TO 260 
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246 PUNCH. I D ( I ) » X ( I ) , Y (I ) » R ( I ) 

PUNCH, ID( I ) »PC( I ) »PI ( I ) * P T ( I ) 
260 I F ( K ( I >-999999.9fl) 250*262 ,26? 
262 CONTINUE 

PAUSE 

GO TO 100 
250 I=T+] 

TF(R< I ) )?00.252»2'>0 
252 RC I l»PIC 1-1 I 

GO TO 200 
C AZIMUTH DISTANCE SUBROUTINE 

10 A»XU+l)-xm 
B«Y< I+ll-YI I > 
PIE«3. 1415927 

11 C«SQR(A*»2+B**2) 
IF(A)15»12.15 

12 IF(B)13»14,14 

13 A2=PIE 
GO TO 18 

14 AZaO.O 
GO TO 18 

15 ANG»ATN(R/A) 

16 IF(A)19*12»17 

17 AZ*(PIE/2.0)-ANG 
GO TO 18 

19 AZ«(1.5»PIE)-ANG 

18 GO TO 1206*100) »ISK4 
FND 



o a 



C RELOCATION OFFSETS / PART 2 

C P. 0. ROBERTS / JULY»1961 

C MIT CIVIL ENGINEERING SYSTEMS LABORATORY 

i ' DIMENSION ID ( 31 ) »Y < 30 ) »X ( 30 ) »R< 30 > ,PC( 30 > »P I ( 30 ) ,PT < 30 ) »OFFD< 3 ) 
C READ TERRAIN , 

K=ID(31)+1 

J=l 

STA=-888888.88 

READ* IBLB»BLSB»DElTB 

IBL=IBLB 
AOS STA=BLSB 

READ. IBLA»BLSA*DELTA 

IF< J-l )409»509*409 
442 BLSB=BLSB+DELTB 

STA=BLSB 

IBL-IBLB 

409 !FUBLA-IBLB>420»410*420 

410 IF(STA-BLSA) 500*500*420 
420 BLSB=BLS A 

IBLB= I BLA 
DFLTB=DELTA 
GO TO 408 

C COMPUTE COORDINATES ON BASELINE 

509 BLSB=BLSB-DELTB 

510 J=J+1 
I = "J 

C GEOMETRY SUBROUTINE 

ISW1=1 

11 



202 ISW2=1 

ISW4=1 
C GO TO A?. PIS 1 

C AZIMUTH DISTANCE SUtiROJTINF 

10 A = X( 1+1 )-X< I ) 
B=Y( I +1 ) -Y ( I ) 
11 C=SOR( A**2+B»«2) 
AZ=0.0 

IF (C-. 0000] )1 8*1 8*1? 
12 ANG=1.0 

TF(B)3»4*4 

3 ANG=-1 .0 

4 B*APSF(A/C) 

P=SOR( ABSF(1.0-B)V2. » 
AZ«SINF(P1 
00 14 L=l*50 
ANGI=P-SINF(AZ) 
IF (ANGI-.1E-05) 15*15*14 
14 AZ=AZ+ANGI 
15 AZ=< ANG*AZ)»2. 
IF(A)17»16*16 

16 AZ=1.5707964-AZ 
GO TO 18 

17 AZ=4.71?3890+AZ 

18 GO T0< 206*605* 660*635) ♦ ISW2 
206 IF( ISV.'l-2)2O7*210*?01 

2C7 AZB=AZ 



ISW1=2 
500 I=J 

PCI=PC(I) 
AZO=A2R 

C CHECK FOR MATCHING BASELINES 

IF( IBL*2-< ID< 1-1 )/100+I0( I)/100))510»515,510 
515 IF(STA-PT( 1-1 ) )442 »520*520 
520 IF ( R ( I )-888 888. 88) 525* 530*530 
525 IF<STA-PCI >530*530»545 

210 ANGI = AZ-AZO 

C ADJUST ANGLES TO LESS THAN 180 DEGRFES 

212 IF(ANGI-3. 1415927) 211*211*214 

211 ANG I = ANGI +6. 2 8 31 8 54 
GO TO 212 

214 ANGI=ANGI-6. 2831854 

IF< I SW4-1) 204*204*690 
C COMPUTE SIGN OF DEFLECTION ANGLE 

204 DS=1.0 

I F( ANGI) 208*2^5* 205 
208 DS=-1.0 

205 ANG=ANGI/2. 

C COMPUTE CURVE TANGENT 

RI=R{I) 

TANG=R I * ( S I N ( ANG ) /COS ( ANG ) )*0S 
AZ=AZ0+3. 141 5927 
C FIND COORDINATES OF PC 

r>IS = TANG 



ISW5=1 
XI = X( I ) 
Y I = Y ( I ) 
C GO TO COORD 1 

C COMPUTE COORDINATES SUBROUTINE 

20 XI=XI+DIS*SIN( AZ ) 
YI=YI+DIS*COS( AZ) 
GO TO (216*217. 540*680) ,ISW5 

216 PCX=XI 
PCY=YI 

C FIND COORDINATES OF CIRCLE CENTER 

AZ=AZO+DS*l . 5707964 

DIS=RI 

ISW5=2 

GO TO 20 
C GO TO COORD 2 

217 XO=XI 
YO=YI 
PCAZ=AZ 

IF(ISW3-2) 550»655*690 
C TANGENT CASE 

530 DIS=PI ( I )-STA 

AZMIT=AZ0+1. 5707964 
ISW5=3 

AZ=AZO+3. 1415927 

XI=X(I) 

YI = Y( I ) 



■ GO TO 20 ^ 2V 



c o c 



© a 



C GO TO COORD 3 

545 I F ( STA-PT { I ) )2 02 ,510,510 
C CURVE CASE 

550 ANG*(STA-PCI)/RI 
X1 = X0 
Yl=YO 

DIS=DS»(-RI) 

AZM I T*AI0+1 . 5707964+DS*A.MG 
AZ=AZMIT 
ISW5*3 
CO TO 20 
C STORE BASELINE COORDINATES 

5A(T'"C0PDX«XI 
CORDY-YI 

C COMPUTE INTERSECTION ON CENTERLINE 

600 I*K 

601 ISW2=2 
60 TO 10 

C GO TO AZ OIS 2 

605 I-I+l 

AZ1-AZMIT 
AZ2«AZ+3 .1*15927 
A«X(I)-COROX 
B«YII»-COROY 
C«SIKf AZ1-AZ2) 
IFIC*. 00001 »632t638»638 
638 0FF«888888.8S 
DlS=-888888.88 



GO TO °610 

C COMPUTE OFFSET FROM BL TO CL AND DISTANCE F ROM' 

632 DO 636 L=l ,2 
DIS=OFF 

0FF=(A*C0S(AZ1)-R*51N< AZ.1 ) )/C 

636 AZ1=AZ2 
610 ANG= AZ 

CLS = Pni)-DIS 
C IS INTERSECTION BEYOND THE PC 

IF(CLS-PC< 1 ) )625. 625,640 
C IS ALIGNMENT FINISHED 

640 IF(R( I )-888888. 88)642. 642,442 

641 1=1+1 

642 IF( ID( I)/100-ID< 1+1)/ 100) 641 ,645,641 
645 AZO=AZ 

K = I-1 

GO TO 601 

C IS ALIGNMENT BEYOND THE PREVIOUS PT 

625 IF(CLS-PT( 1-1) )628, 626,626 

626 IF( I-K-l 1698.698.697 
697 K=K+l 

GO TO 698 
628 IFU-K-1 J650. 442.650 
C INTERSECTION ON CURVE 

650 1=1-1 
1SW3=2 
ISW4=1 
GO TO ?10 



C GO TO LAST PART OF GEOM 2 

655 A*XO-CORDX 
B«YO-CORDY 
ISW2-3 
GO TO 11 

C GO TO LAST PART OF AZ D1S 3 

C COMPUTE CURVE OFFSET 

660 ANO-AZ-AZMIT 

A«C»COS(ANG> 

B-C*SIN( ANG) 

OFF»A-DS»SGR C R 1**2-B**? ) 

C GET COORDINATES OF INTERSECTION POINT 

XI=CORDX 

Yl«CORDY 

AZ-AZMIT 

DIS«OFF 

ISW5»4 

GO TO 20 
C GO TO COORD 4 

C GET AZIMUTH OF RADIUS 

680 A-XO-XI 

8-YO-YI 

ISW2-4 

GO TO 11 
C GO TO AZ DIS 4 

C COMPUTE ANGLE SUfiTENDED 

685 ANGI»AZ-PCAZ 

ISW4=3 



GO TO 212 
C GO TO GFOV 3 

C COMPUTE CLS AND LOCAL AZIMUTH FOR CURVE iNTLKStCT ION 

690 CLS = PC ( I ) +ANG I *R I *0S 
ANG=AZO+ANGI 

698 ANG =SIN(AZMIT-ANG) 
IF(SENSE SWITCH 1)699*700 

699 PRINT* IBL»STA.CLS»OFF»ANG 
GO TO 442 

700 PUNCH* IBL*STA»CLS*0FF »ANG 
GO TO 442 

END 



o o 



o 



REFERENCES 



(1) "Using New Methods In Highway Location", Paul Roberts, 
Photogrammetrtc Engineering, June, 1957 

(2) Digital Terrain Model Approach To Highway Earthwork Analysis, 
C. L. Killer, MDPW - BPR Research Project Report, Massachusetts 
Institute of Technology, August, 1957 

(3) The Skew System For Highway Earthwork Analysis, C. L. Miller 
Massachusetts Institute of Technology, MDPW - BPR Research 
Project Report, September, 1957 

(k) Earthwork Data Procurement By Photogrammetrlc Methods, 

C. L. Miller, and T. H. Kaalstad, Highway Research Board 
Bulletin 199 

(5) "The Digital Terrain Model - Theory And Application'*, 

C. L. Miller, ami R. A. Laflamme, Photogrammetrlc Engineering, 
June. 1958 

(6) Digital Terrain Mode T System Manual, C. L. Miller, and 
R. A. Laflamme, MDPW - BPR Research Project Report, 
Massachusetts Institute of Technology, December 1958 

(?) The Digital Terrain Model System of Providing Highway 

Location And Design Information, C. L. Milter, Proceedings 
Of The American Association Of State Highway Officials, 
Committee Meeting, San Francisco, December 1958 

(8) Preliminary Report on The Digital Terrain Data Recorder, 

C. L. Miller, and E. P. Gladding, MDPW - BPR Research Project 
Report, May 1959 

(9) Zone Cost Evaluation Program EA-2, C. L. Miller, and 

D. E. Weisberg, MDPW - BPR Research Project Report 
Massachusetts Institute of Technology, October i960 

(10) Digital Terrain Model System DTM II Manual, P. 0. Roberts, 
C. L. Miller and R. A. Laflaimie, MDPW - BPR Research Project 
Report, Massachusetts Institute of Technology, November I960 

(11) A New Technique For The Prediction Of Vehicle Operating Cost 

In Connection With Highway Design, A. S. Lang, and D. H. Robbins, 
MDPW - BPR Research Project Report, Massachusetts Institute Of 
Technol ogy , Nov ember 1 960 

(12) Use Of Digital Computers In Land Cost Evaluation, D. E. Weisberg. 
Unpublished S. M. Thesis, Massachusetts Institute Of Technology, 
June 1961 

(13) Highway Location Evaluation, P. 0. Roberts, Proceedings Of The 
Electronics Committee Of The American Association Of State 
Highway Officials, Boston, August I96I 











C kLLoCATI Gim OFFSETS / PAR I . 1 

C P. KJdLKTo / JULY » 1961 

C MT CIVIL LN OlUccRIN C SY oTEiMS LAbOKA.TO.HV 



c 


w>w 1 UN FOR TYPED OUTPUT / MAR 14»95 / TABS 22*34 


46 58» 


70»82 




Ulr.tKoJUN ID I 31 ) » » I 3U ) » A I 3U 1 »K I 30 ) tPC ( 30 J iP I ( 30 ) 


tPT ( 3Q1 


iOFFD(3) 


■ ■ - - - 


DI.mc.NS ION D( 30 ) » T ( 30 ) »|£( 30 ) 






c 


T M T T T A 1 1 7 U" 






100 


I =0 






c 


READ ROUTINE 








D- 262 N=l»2 








I D ( 3 1 ) = J 






• 


DO 1Z<J J=l»30 









.1=1+1... . ... 








READ» ID(I)»X(I)»Y(I)»R(I)*D(I>»T(1)»E(I) 














120 


C0NT I NUE 






140 


I = I.-J+1 






200 


I SW 1=1 
I S »'» 3 = 1 








RI =R { I ) 








P I ( I ) =R I 








PC ( I ) =RI 








PT ( I ) =R I 








C = • 








T(I)=0. 




/ 

j 


5 '1 1 


I S 1 

uU IU I U 






r c 








248 


i O W J. — £. 








A M C 1 ( 1 — I i n 








MLHI\Lj— U • U 








cn Tn 9^9 






5 1 n 








91 n 
c. 1 e. 

211 










cn T n 9 i 9 






214 


ANG I = ANG I— 6.2831854 














2^8 


Q3 = , ~l #UuG0GQG 








r.n Tn 9 l i ^ 






2u9 


D S = 1 • w U u 






2^5 








C 


FIND w.HICH IS GIVEN 








T F ( k ( T 119 5 1 «91£».991 






216 


!C( D( T 11951 .517.910 






Cxi 








7 1ft 


IFIPI1U551.555.551 






C 








219 


R( I )=5729.5780/D( I ) 








60 10 221 






C 


TIT T TC r t i/PM 

IF 1 IS GIVEN 








K \ I ) - 1 l I J *tUo I ANG ) / o I N I ANG ) *D«s 
GO TO 221 . 






C 


IF fc IS GIVEN 






223 


R ( I ) = E ( 1 ) / ( 1 • /COS ( ANG ) ~1 • J 








GO TO 221 








IF NONE ARE GIVEN 






222 


D(I)=0.0 

GO TO 215 






221 


D ( I )=t>729.5760/R< I ) 






215 


T( I )=R( I )*(SIN(ANG)/COS(ANG) )'*DS 



















O i- 1 i ; f~ A T 1 / ' t < C c c_ t_' T r / n «-i t -i 













C 

PAC ^ 


P. 0. KuoLk lb / JULY»lv61 "'" ' "~ ~ 






PC ( I ) =PI ( I )-T ( I ) 

ALEhG=ANGI*R( I )*DS 








OI.il-.MolON ID (31 ) »Y < 30) ,X (30) 30) ,PC< 30) »PI < 30 ) »PT ( 30 ) .0FFD ( 3 ) 

K- J - A i i T 1- 13 P A T M 






PTC I )=PC (I J+ALENG 

Ei I ) = ( l./C0S(ANG)-l. )*R{ I ) 








K=ID( 31 )+l 

<J — x 




7 "\ P 
fc. 3 C. 


ANG I 0=ANGI*57. 295780 
AZ0UT= AZ*57. 295780 








iTA*-8aUb8C.U8 

f\FAr)»Ii-l ^^rtl 'iA.I^-l Tli 
t\^MU9 li_iL.09DU Oj3 > U U U 1 y 






IF(i>ENSE iftlTCH 1)233»234 
PRINT » 10 ( I ) »X( I ) ,Y( 1 ) ,R( I ) »D ( I ) 


»T( 1 )',!£( I ) 






IUL=Ii3Ld 






PRluT»It>( I ) tPCU ) »PI ( I ) »PT (I ) »AZ0UT »ANGI0»ALENG 
GO TO 235 






KEAJ » 1 uLA* bLSA »DELTA 

r c ( 1—1 »/i^o »i\iio ./wif) 




2 34 


PUNCH* ID ( I ) »X( I),Y{ I ) ,R( I ) ,D( I ) »T( I ) ,F( I ) 
PUNCH* IOU) »PC(I) »PI ( I ) »PT til »AZOUI ,ANGIQ»ALENG 






o 1 A-dLOu 




23i> 


PI I 1*1 isPT t I \ — T t t \ -k.r 

AZC=AZ 






4w9 


I BL= I ULb 

ir l 10LA~luLuf fcU>f iUtfzU 






i * * ~ i 

I.FtRC I )-888888.87) 242*244*244 






41b 
420 


IF(iTA-BLSA)500,500»420 

OLOu-*DLOM 




242 


fdhl f T i\UJF 

GO TO 201 








IBLB=I3LA 

U C L * o — U u L, 1 P\ 




244 


PCI 1 1*PI < I > 
PTMl-PHU 








GO TO 408 

rOiViPLiTf 1 ' rOORfi T WAT F t, f,M ujtfi IMC 




4.HD 


IF(SEN£E SWITCH 1)245*246 
PRUmT. 1D( I ) »X( I ) »Y< I ) »R< I > 


— _ 




5u9 
510 


BLSB=BLSB-DELTB 
J=J+1 






PRINT* IO(ir»PC(I)»PI(I)*PT(I) 
GO TO 260 








I=J 

\J W ' 1 u> 1 l\ 1 o uy l\ w W | i lit 




246 


PUNCH* 10(1 ) »Xt I ) »Y ( I ) »RU ) 
PUNCH* ID( I > »PC( 1 ) »PI U )»PT(I) 


: ~ 






ISW1=1 
1 = 1—1 




26v> 


IF(Rt I )-999999«9a)250»262»262 
CONTINUE 






2U2 


ISW2=1 • ~ 
ISV.'3 = 1 






PAUiE 

GO TO 100 






Q 


ISW4=1 

(~if) Tn A7 ni<; i 




250 


.1*1 + 1, 

IFCRU>)200*252»200 






C 


AZIMUTH DISTANCE SUBROUTINE 

A = X ( I + 1 \ —X < I 1 




252 


R 1 1 )*P1 1 I— I 1 

GO TO 200 






11 


B=Y( I+l)-Y< I ) 
C=S0R ( A**2+B**2 ) 




C 

10 


AZIMUTH DTSTAMrp ^IJRRfHiT I NF 
A=X( I + 1 ) -X ( I ) 








AZ=u.O 

I F ( C-. 00001 ) 18 » 18 » 12 






8*YM + 1)-Y(I) 
PIE=3. 1415927 






12 


ANG=1.0 

I F ( Li ) 3 »4 »4 




11 


C»SQRIA**2+B**2> 
IF(A)15.12.15 






3 


ANG=-1.0 

B=AbSF ( A/C ) 




12 
13 


IF(B)13»14»14 
AZ*PIE 








P=SQR CABSF ( 1 ,0-B ) /2 • ) 

AZsS I NF (Pi 

n*. — ^ 1 111 \ i / 




1 A. 


GO TO 18 
AZ»0.0 








DO 14 L=l»50 

AM^ ! -D — ^ IMP / A 7 1 
MINO 1 — r—Ca 1 INr I At J 




1 5 


GO TO 18 
ANG=ATN( B/A) 






14 


IF(ANGI-.1E-05)15»15 ,14 
AZ=AZ+ANGI 




16 
i t 
l / 


IF(A) 19*12*17 
AZ=(PIE/2.0)-ANG 






15 


AZ=(ANG*AZ)*2. 

IrlMJJ. / » lo » io 




19 


GO TO 18 

AZ*(1.5»PIE)-ANG 






16 


AZ=1.5707964-AZ 

bU 1 U 10 




18 


GO TO (206*100 >»ISW4 

END 






17 
18 


AZ=4.7123890+AZ 

GO T0< 2o6»605,660,685 ) »ISa2 












206 
2w7 


IF( IStvl-2)207,210,201 

AZB=AZ 














1 = 1 + 1 <IQ 

ISW1=2 J 












500 


PCI-PC(I) * 






o - o ■ 






G ■ 





© a 



o © 



■S AT CUI NG B AS E L I N E S_ _ 



515_ 

52u 

525 



21o 

Tl2 _ 
211 



214 



AZ0=AZb 

CHE CK FO R 

IFUBL*2-( IDC I -1 ) /100 + l5( I )7lOO ) ) 310,515,510 

Jfl^TA-PTiI^JJ.M_2_».5ZQ>i)2y^ 

IF(R (I)- 88 8868. 66)525,530,530 

IFt QTA- PCI )530,53 0,545, 

ANGI = AZ-AZO 

ADJUST ANG LES TO LESS THAJi_180 DEGREES 
IF(aNGI-3. 1415927) 2 i 1,2 i 1 ,2 14 "~ " 

ANGI=ANGI+6.2831654 

GO TO 212 

ANGI=ANGI-6.2831854 



PAT. 



DIS=Do*(-RI ) 

AZI«;i T = AZO+1.57_07964 + DS*ANG _ 

az=az;-;it 

Ioi-.5 = 3 _ 

GO TO 20 

I oTOKL BA SELINE COOKD lNAl to 

540 COKDX=Xl 

_ ccrdy=yi 

i' "~ compute ~ intersect ion on centerl i ne 

_6O0 I_=K 

6ul ISU2=2 

GO TO 10 



c 


I F I ISsa.4-1 ) 204 ♦ 204 »690 

COMPUTE SIGN OF tiEFLF CT T OK ANCIF 


O \j \> 


co Tn a 7 nit, o 

1-1+1 




DS=1.0 

IF (ANGI J 2U8»205»2fl5 




A J 1 — A 7 iV. I T 


208 


DS=-1.C 
ANG=ANGI/2. 


: 


a=x ( i ) -f osnx 

O-T I I J — V-UKl> T 


c 


COMPUTE CURVE TANGENT 
RI=R ( I 1 








TANG=R1* ( S IN ( ANG ) /COS ( AiiG ) )*DS 
AZ=AZ0+3. 1415927 


638 


flFFsHHHHRft- Rtt 
urr-oooouOi 00 

J10--OOOBt)o.Ofl 


C 


FIND COORDINATES OF PC 




CO TO AT n 




DIS=TANG 




COMPUTF QFF^FT FROM HI TO Tl ANiM flKTiwrr FI?r>M Dl 




ISW5=1 
XI=X( I ) 


w J C 


DIS=0FF 


c 


YI = Y ( I ) 

GO TO COORD 1 


636 


AZ1=AZ2 


c 


COMPUTE COORDINATES SUBROUTINE 


610 


ANG= AZ 


20 


XI =X I+DI S*S IN( AZ ) 








Yl=YI+DIS*COS(AZ) 


C 


IS INTERSECTION BEYOND THE PC 




GO TC ( 2 16 » 2 17 » 540 » 6 80 ) » I SW5 




I F ( C" I S— PC ( I 1 lA?S.A?h.AAi'i 


216 


PCX=X I 


Q 






PCY=YI 


64u 




c 


FIND COORDINATES OF CIRCLE CENTER 


641 


1 = 1 + 1 




AZ=AZO+DS*1.5707964 


642 


IF( ID( I)/100— ID( I + l)/100)641f645» 641 




DIS=RI 


645 


AZ0=AZ 




ISW5=2 




K=I-1 




GO TO 20 




GO TO 601 


c 


GO TO COORD 2 


C 


Io AL I GKi'iENT BEYOND THE PKEVIOuS PT 


217 


X0=XI 


625 


IF ICLJ-PT( 1-1) ) 628 ,626,626 




Y0=Y I 


626 


IF( I-K-l )698»698,697 




PCAZ=AZ 


697 


K = K+1 




IF( ISrt3-2) 550,655»690 




GO TO 698 


C 


TANGENT CASE 


628 


IF( I-iv.-l )650»442,650 


530 


DIS=PI ( I )-STA 


C 


I NT tRSECT I ON ON CURVE 




AZM1T=AZ0+1. 5707964 


650 


1 = 1-1 




ISw5=3 




ISW3=2 




AZ=AZ0+3. 1415927 




ISW4=1 




XI=X(I) 




GO TO 210 




YI=Y< I ) 


C 


GO TO LAST PART OF GEOM 2 




GO TO 20 


655 


A=XO-CORDX 


c 

545 


GO TO COORD 3 

I F ( oTA-PT ( I ) )202p510»510 


V 


B=YU-C0RDY 

ISW2=3 *P 


C 


CURVE CASE 




GO TO 11 


550 


ANG=(STA-PCI ) /RI 


C 


GO TO LAST PART OF AZ DIS 3 




XI = XO 


c 


COMPUTE CURVt OFFSET 




YI=Y0 


660 


ANG= AZ-AZM I T 



A=C*CUS(ANG» 

B=C»SIN( ANG ) 

QFF=A-D_*S_R(RI**2-U**2 ) 

Q_T COORD I NATE S OF I MI£R^LCLlC)N_£.ai^.L 

XI=CORDX 

YI=COKDY 

A_*AZ>iIT 

_IS=OFF 

' I S»*5*4 

GO TO 20 _ 

C GO TO COORD 4 

C GET AZIMUTH OF RADIUS 

680 A*XO-XI 

3-YO-YI 

ISK2=4 

GO TO 11 

C GO TO AZ DIS 4 

C CQrtPUTE. ANGLE SUBTENDED , 

685 ANGI=AZ-PCAZ 

lSrt4«3 

GO TO 212 

C GQ TO GE QM 3 

C COMPUTE CLS AND LOCAL AZIMUTH FOR CURVE INTERSECTION 

690 CLS»PC( I )+ANGI«RI*D. . 

ANG*AZO+ANGI 

698: ANG IN ( AZMI T-ANGl 

IFUEWiE -WITCH 1)699.700 
699 PRlNT»I.L»-TA»CL,»OFF»ANG 

GO TO 442 

7U0 PUNCH»l-Lf-TA>CLS.OFF»AN- 

GO TO 442 

; , H : £Np, , . : ■ , , 



o 



c 



O 



